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We consider the influence of a resonant control field on weak localization of light in ultracold atomic ensem- 
bles. Both steady-state and pulsed light excitation are considered. We show that the presence of a control 
field essentially changes the type of interference effects which occur under conditions of multiple scattering. 
For example, for some scattering polarization channels the presence of a control field can cause destructive 
interference through which the enhancement factor, normally considered to be greater than one, becomes less 
than one. 



1. Introduction 

Interference effects under multiple light scattering in optically thick disordered media have been 
subject to intense investigation for almost three decades. Coherent backscattering (CBS), which 
is closely related to weak localization, is one of the striking examples of the types of effects 
which can occur. CBS manifest itself in an enhancement in the intensity of light scattered in the 
nearly backwards direction. Wave scattering in this direction along reciprocal, or time-reversed, 
multiple scattering paths preserves the relative phase, which results in constructive interference. 
The first detailed observations and analysis of this effect for light were made in (I3)-(Q) and by 
now CBS in solids and liquids has been investigated in detail (Q)-(0)- 

Observation of CBS of light in atomic gases is complicated because atomic motion causes 
random phase shifts of scattering waves which are different for reciprocal paths. For this reason, 
the weak localization has been observed only for a cold atomic ensemble. The first experiments 
on CBS in ultracold atomic samples prepared in a magneto-optical trap shown that 

there are a range of interesting features of this phenomena compared with CBS observed earlier 
for solids and liquids These features, which have their origin in the atomic nature of 

the scatterers, could not be understood and quantitatively described by approaches developed 
previously for classical scatterers such as powders and suspensions. 

Detailed treatments of CBS, takin g in to account the quantum nature of atomic scatterers, 
was developed by several groups ©-Hi)- In these papers different aspects of weak localization 
were considered. Particularly it was shown that one can have a strong influence on the observed 
interferences by reinforcing those scattering channels which lead to interference and suppressing 
those which do not. Such desirable effects can be realized for example by polarization of atomic 
angular moments. By means of optical orientation effects it is possible to collect all atoms in one 
Zeeman sublevels ensuring the fulfillment of optimal interference. Similar effects can be achieved 
by applying an auxiliary static magnetic field and tuning the frequency of the light in such a 



way that light would interact only with the desire Zeeman sublevel (Lo). 

Experiment {If} has shown that a static external field influences the process of multiple 
scattering and the associated interferences. At the same time it is well known that control of 
optical properties of matter by means of auxiliary quasiresonant electromagnetic fields is much 

ISSN: 0950-0340 print /ISSN 1362-3044 online 

© 200x Taylor & Francis 

DDI: 10.1080/0950034YYXXXXXXXX 

http:/ /www.informaworld.com 



January 13, 2013 14:52 Journal of Modern Optics CBSvsEITSubmitted 
2 

more effective than by a static one. Electromagnetically induced atomic coherence changes the 
optical properties of atomic samples in sometimes dramatic ways, and is responsible for such 
effects as population trapping, electromagnetically induced transparency (EIT), "slow light" , 
"stopped light," to name a few Hi)-®. 

In this paper we are going to analyze how a resonant control field influences coherent backscat- 
tering of light. In particular, we consider CBS effects under conditions of electro-magnetically 
induced transparency. Our attention is focused on the spectral dependence of the relative am- 
plitude of the backscattering cone in a steady-state regime, i.e. on the dependence of dthe 
enhancement factor on frequency of scattered light which assumed to be monochromatic. We 
show that a control field can essentially change the type of interference and even can cause de- 
structive interference. For some scattering polarization channels and for some detunings of the 
probe light, the enhancement factor can be less than one. With the developed knowledge of the 
spectrum we also consider the dynamics of CBS in the case of pulsed probe radiation. 

2. Basic assumptions 

One of main quantitative characteristics of CBS is the enhancement factor, which determines the 
relative contribution of interference effects to the total scattered light intensity. In experiment it 
is measured as the relative amount of light intensity scattered into a given direction inside CBS 
cone to the background intensity registered outside the cone. In theory it is more convenient to 
evaluate the differential cross section of scattering from the input to outgoing mode and calculate 
the enhancement factor as the ratio of this cross section to its non- interfering part. 

Our theoretical approach allowing calculation of the differential cross section of light scattered 
from optically thick ultracold atomic ensembles is described in details in a series of papers (fiil ). 
(fi^ )- (p^ ). This approach is based on a diagrammatic technique for nonequilibrium systems and 
allows us to obtain separately both interfering and noninterfering parts as a series over a number 
of incoherent scatterings. 

The generalization of this approach to the case of the presence of a coherent control field was 
made in (0) and [2^). In Appendix A we show, as an example, the double scattering contri- 
butions to the differential cross section. On the basis of these and similar expressions for higher 
order scattering contributions, we calculate here the spectral dependence of the enhancement 
factor. We consider probe light scattering from ultracold clouds of ^^Rb atoms, prepared in 
a magneto-optical trap after the trapping and repumping lasers and the quadrupole magnetic 
field are switched off. All atoms populate the F=l hyperfine sublevel of the ground state, while 
the distribution over Zeeman sublevels is uniform. The spatial distribution of atoms is assumed 
to be spherically symmetric and Gaussian. For the typical conditions of the trap, the Doppler 
width is many times smaller than the natural line width of the excited state and the interatomic 
distances on average are much larger than the optical wavelength (dilute medium). This allows 
us to neglect all effects associated with atomic motion, and atomic collisions. 

Probe radiation is quasi resonant with the F = 1 ^ F' = 1 transition of the Di line (see Fig. 
1) and its polarization can be arbitrary. However, for definiteness we will consider right or left- 
handed circularly polarized light. This light is assumed to be weak; all nonlinear effects connected 
with the probe radiation will be neglected. In our calculations this field will be taken into account 
only in the first non-vanishing order. Besides the probe light, the atomic ensemble interacts with a 
coupling, control field. In this paper we will consider this field tuned to exact resonance with F = 
2 — )• F' = 1 transition. Its amplitude is determined by the Rabi frequency Qc = 2 1 Kim' I) Kim' are 
the transition matrix elements for the coupling mode between states In) and \m') = \F, m'), which 
we define with respect to the \m') = \F = 2,m' = —1 > — t- |n) = \F' = l,n = > hyperfine 
Zeeman transition. Other transition matrix elements are proportional to Qc and algebraic factors 
depending on corresponding Clebsch-Gordon coefficients. The polarization of the control field is 
assumed to be right handed and circular. The detuning of the probe laser frequencies from the 
corresponding atomic resonances is assumed to be much less than the hyperfine splitting. This 
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Figure 1. Excitation scheme for observation of the EIT effect in the system of hyperfine and Zeeman sublevels of the 
Dl-line of ^'^Rb. The coupling field is applied with right-handed circular polarization to the _F = 2 — ^ F' = 1 transition. 
The probe mode is applied to _F = 1 — > _F' = 1 transition and can cause different excitations depending on its polarization 
and propagation direction. 



circumstance, along with the relatively large hyperfine splitting in the excited state, allows us 
take into account only one hyperfine sublevel F' = 1 oi this state. 

In this paper we will not focus our attention on the angular distribution and shape of the 
CBS cone, but will instead restrict ourselves by consideration of exact backscattering only. The 
ground state hyperfine splitting of ^'^Rb is about 6.83 GHz, so Rayleigh and elastic Raman 
scattering (associated with the F = 1^F' = 1^F = 1 transition) is well spectrally separated 
from inelastically Raman scattered light {F = 1 ^ F' = 1 ^ F = 2). We will assume spectral 
selection under photo detection. The inelastic Raman component is assumed to be not registered. 
Depending on the type of polarization analyzer used for photodetection, four main polarization 
schemes can be considered — )■ H^, — )• — )• H^, — )• Here H± represents 
the helicity of the input and outgoing light. Note that, despite homogeneous population of the 
Zeeman sublevels, the susceptibility tensor becomes essentially anisotropic due to the presence 
of the coupling field ( In such a case the enhancement factor depends not only on the 
relative polarizations of the input and output light but also on their specific types. 



3. Results 

The results of calculations of the spectrum of the enhancement factor for a weak control field 
are shown in Fig. 2a. The calculations were performed for an atomic cloud with a Gaussian 
radius equal to = 0.5 cm; maximal density is no = 3.2 • lQ^^cm~^. The Rabi frequency of 
the control field is Jlc = O.SF. For a weak control field we observe a not particularly surprising 
behavior of the spectrum. Against a background typical for the spectral dependence of CBS 
effect we see narrow spectral gap which has its origin in the decreasing of the optical depth of 
the cloud caused by the EIT phenomena. Under the EIT effect the probability of higher order 
scattering relatively decreases compared with single scattering and all interference effects which 
connect with multiple scattering are reduced. The width of the gap in the spectrum is about the 
width of the transparency window determined by the EIT phenomenon. We point out only that 
there is a certain difference in the width for different polarization channels caused by the above- 
mentioned optical anisotropy of the atomic ensemble. The situation changes dramatically when 
the control Rabi frequency becomes comparable to or larger than the spontaneous decay rate. 
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In Fig. 2b, which is calculated for J7c = SF, in addition to the more noticeable anisotropy, an 
essential transformation of the spectrum takes place. The range of the structure of this spectrum 
connects with the difference in the Autler-Townes splitting for different Zeeman transitions. The 
latter is caused by different dipole moments of the corresponding transitions and consequently 
with different Rabi frequencies for them. The maximal value of the enhancement factor for 
polarization channels with changing helicity is almost the same as for weak control field but 
for the case with preserving helicity we see qualitative modifications. The main one of these 
is that, instead of constructive interference for some spectral regions, we observe destructive 
interference. In these regions, for channels -ff_ — t- and — t- the enhancement factor 
becomes less than unity. That is, in place of a CBS cone we have a CBS gap, or anticone. In 
spite of the relatively small value of the gap it seems physically important because CBS or 
weak localization itself in its "traditional" interpretation connects with time-reversal and always 
causes enhancement in back scattering. 




A similar effect for scattering of polarized electrons in a solid was shown in (j^j ). Following 
this, antilocalization in electron transport has been widely studied in a range of physical systems 
where electrons interact directly with magnetic impurities and where the spin-orbit interaction 
is important (24, [H). The possibility to observe destructive interference in the case of light 
scattering from ultracold atomic ensemble was predicted for the first time in (j7^. This effect 
was explained in (f73 ) by the hyperfine interaction in atoms and by possible interference of 
transitions through different hyperfine sublevels of excited atomic states. In the case considered 
now, the mechanism of antilocalization is different. This is emphasized by Fig. 2b, in which 
the results are obtained for only one excited state sublevel F' = 1 without taking into account 
possible hyperfine interactions. 

The fundamental possibility of destructive interference under multiple scattering connected 
with the coupling field is illustrated in Fig. 3. Here we show double backscattering of a positive 
helicity incoming photon from the probe light beam on a system consisting of two ^''Rb atoms; 
the exit channel consists of detection of light also of positive helicity. Here the double scatter- 
ing is a combination of the Raleigh-type and Raman-type transitions. In the direct path the 
scattering consists of a sequence of Rayleigh-type scattering in the first step and of Raman-type 
scattering in the second one. In the reciprocal path, Raman-type scattering occurs first, and 
the positive helicity photon undergoes Rayleigh-type scattering in the second step. There is an 
important difference in the transition amplitudes associated with Raleigh process for these two 
interfering channels. Indeed, in the direct and reciprocal path the scattered mode is coupled with 
different Zeeman transitions. In the absence of a coupling field these transitions have the same 
amplitude and we always have constructive interference. The coupling field essentially modifies 
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the scattering process and this modification is different for difi'erent Zeeman transitions. For def- 
inite frequencies of probe hght, the scattering amphtudes connecting the direct and reciprocal 
scattering channels can be comparable in absolute value but can have phases shifted by an angle 
close to TT. In this case these two channels suppress each other. The considered example is not 
the only one, and appears only in the double scattering channel. There are some others which 
lead to constructive interference. In the higher scattering order the situation is similar. That is 
why in our results we have only partial destructive interference and for the considered param- 
eters only a small suppression of backscattering. Note that the antilocalization effect discussed 




Figure 3. Diagram explaining the antilocalization phenomenon in the example of double scattering in the helicity preserving 
scattering channel. Only the interfering transitions initiated by the probe are shown and the presence of the control mode 
is specified in Fig. 1. 

in (idi ) takes place for essentially nonresonant light and consequently is difficult for observation, 
because it is challenging to prepare clouds with large optical depth for such radiation. The effect 
considered here takes place for radiation which is not far from resonance. It makes this case 
more promising for experimental verification. Note also that similar antilocalization phenomena 
can be observed in the case when the control field is absent but the probe radiation is strong 
enough to make the Autler-Townes splitting noticeable (H^). 

Consider now how these peculiarities in the spectrum manifest themselves in time dependent 
CBS in the case of a pulsed probe light. Here the most interesting effects are observed for a 
weak control field when peculiarities in the spectrum are sharp. Different spectral components 
of the input pulse scatter differently and it is the reason for the essential transformation of the 
spectrum and consequently the time profile of the pulse. Scattered light has a deficit of near 
resonant photons for which the EIT mechanism works best. The gap in the spectrum causes light 
beating effects which are different for scattering of different orders. In Fig. 4 we show shapes 
of pulses which undergo single (4a) and double (4b) scattering for the if+ — t- polarization 
channel. These graphs are calculated for a cloud with no = 3.2 • 10^^cm~^ and tq = 0.5 cm. The 
input pulse has a Gaussian shape / = /q exp(— t^/r^), the length of the pulse is T^200r^^. 

The shapes of single and doubly scattered pulses are essentially different. Most important 
is that the intensity peaks for them are separated in time. When single scattering intensity is 
maximal double scattering is very small and vice versa. This effect arises through the relative 
time scales associated with single scattering in comparison with double scattering, as illustrated 
in Fig. 5. For comparison the input pulse shape is also shown. 
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Figure 5. Time dependence of tlie enhancement factor in the case of pulsed probe radiation. Input pulse shape (in arbitrary 
intensity units) is shown in gray 
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Appendix A. The cross section of the scattering process 

The scattering process is typically described by the differential cross section for light scattered 
into the solid angle 17. If the system is dilute but optically dense than the total outcome can be 
expressed as the sum of the contributions of the subsequent scattering events, see 

da ■^-y da ■^-^ daat v~~^ daatc 

dn^ ^^dn^ ^ In ^ dn ^ ' ' ' ^ ^' 

a=l a^b aji-h^c 



Through macroscopic electrodynamics based self-consistent averaging in the dilute disordered 
system, this series is converging rapidly such that any events of possible recurrent scattering 
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(when the atomic numerator gives a = c etc.) can be safely ignored. We consider as an example 
the calculation scheme in the second order which can be straightforwardly generalized on the 
higher orders. For the specific scattering channel in the backward direction the second order 
contribution from two randomly selected atoms is given by the sum of two terms 

dn dn dn ^ ' 

where the first term is normally called a " ladder" term describing the contribution of successive 
double scattering of a photon of frequency uj, wave vector k, and polarization e to the outgoing 
mode a;',k',e' on atom a = 1 first and atom 6 = 2 second. A similar term contributes to 
the expansion ()Aip for the reciprocal path with a = 2 and 6 = 1. The second "crossed" term is 
associated with interference between these scattering events. For the sake of clarity let us specify 
the polarization of the incoming photon by Cartesian unit vector with ^ = x,y (k || z) and 
of the outgoing photon by with v = x' ^y' (k' || z'). 

Then scattering in the particular polarization channel from to ej, is described by successive 
double scattering in the "ladder" contribution and is given by 

^'^12^ \ ^ \ ^ U}\2 ^^'^ 1 



E E 



dn ^ ^ c4 c4 rl 



ij^^yl i,j,k,l 



12 



xXj,i(oo,r2,a;')a.fc' '(a;i2-ki2V2)XH(r2,ri,w^jQ;;^.^ ' (a;-kvi)Xj^(ri, -oo,w) 

X^j{oo, r2, uj')a?^ (wi2 -ki2V2)X^j(r2, ri, ^12)0^^' ™' (w-kvi)Xj^(ri, -oo,w)J (A3) 

where the photon's output and virtual intermediate frequency are affected by both the Raman 
and residual Doppler effects 

W12 = W - UJm'^rm + (ki2 - k)vi 

= UJ12 - ^m'im^ + (k' - ki2)v2 (A4) 

where the following notation is used: vi and V2 are the velocities of the first and the second atoms; 
ki2 and k' are the photon's wave vectors in the intermediate and outgoing modes; ryi = |r2 — ri| 
is the relative distance between the atoms. For the intermediate mode the wave vector is given 
by ki2 = (u;i2/c)(r2 - ri)/|r2 - ri|. 

There are several important ingredients contributing to the cross section. The first is the 
scattering tensor, which is given by 



a 



(m" m) , 1 {dl)m"n{dj)nTn ( ^ \^nm'\ 1 I /' A 

^ A + ir/2 \ ~ A + iT/2 Ae- A + S„„,(A) f ^ ' 

m'{n),n ^ / . 



With reference to Fig. 1 the self-energy radiation term S„m'(A) = |K„rrt'P/ (A-|-zr/2) is ex- 
pressed by transition matrix elements Vnm' for the control mode coupling the non-populated 
states |n) and \m') = \F,m') the frequency detunings Ac and A are respectively the offsets of 
the probe uj and control Wc from the resonance A = w — ojf'f {F = 1) and Ac = Wc — ojf'F 
{F = 2); r is the natural radiation rate of the excited states; {dj)nm represents the dipole tran- 
sition moments between the lower \m) = \F,m) and upper |n) = \F',n) states. The scattering 
tensor determines the amplitude of the single scattering event for either the elastic or inelastic 
channel accompanied by the atomic transition from the state \m) = \F,m) with F = 1 to the 
state \m") = \F,m") with F = 1 (Rayleigh channel) or F = 2 (inelastic Raman channel), see 
Fig. 1. 
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Before, between, and after two successive scattering events the photon propagates in the bulk 
medium in accordance with the macroscopic Maxwell theory. The transformation of its amplitude 



is relevantly described by the formalism of the Green's propagation function, see ilA . \2j\ . \2£). 
The matrix functions Xfc;(r2, ri; w) perform the slowly varying amplitudes associated with the 
photon retarded-type Green's function. They follow the transformation of the light amplitude 
and polarization between points ri and r2. One of these points can approach infinity for either the 
input or output scattering channels. The subscript tensor indices k and / are always associated 
with the specific frame where the propagation direction r2 — ri is associated with z-axis. The 
components k and I belong to the orthogonal plane z and we enumerate them by 1 (x-axis) and 
2 (y-axis). In this basis the Green's function slowly varying components can be expressed as 
follows 

Xn(r2, ri; = e'^'^''^'''^) (cos <^(r2, ri) - sin 0(r2, ri)n,) 
X22{r2,ri;uj) = e^'^°(r„rO ^^^^ 0(r2, n) + sin 0(r2, ri)n^) 
Xi2(r2, ri; tj) = ie'^oi^^'""') sin0(r2, ri) {uy + in,) 

X2i(r2,ri;^^) =ie*'^«('-^''-i)sin0(r2,ri)(nj,-m,) (A6) 

where 

0o(r2,ri) = / xoir,u))ds; 0(r2, ri) = / x(r,w)(is (A7) 

^ Jri ^ Jri 

performs the phase integrals along the path s linking the points ri and r2 such that r = r(s) in 
the integrand. The parameters of these integrals are expressed by components of the dielectric 
susceptibility tensor of the medium. 

Because of the specific symmetry of the problem, see Fig. 1 it is convenient to define these 
integrals in the alternative basis set of circular polarizations. It can be defined by the following 
expansion eg = e^, e±i = =p(e^ibiej/)/\/2 in the Cartesian frame and it requires co/contravariant 
notation in writing the tensors products, see ^2^. The spectrally dependent parameters x(r,<^), 
Xo{'^) ^-ncl symbolic vector ^(w) perform the expansion coefficients of the susceptibility tensor 
projected on the plane orthogonal to r2— ri in the basis set of the Pauli matrices ~^ = (cr^, (Ty,az)- 
This tensor projection thus can be written 

x/(r,a;) = Xo(r,a;)V' + (^(r,^) • (A8) 

where q,q' = ±1. To find the parameters of the phase integrals the components Xg'^ iii the left 
side of the expansion (jASP should be explicitly expressed by the components of the susceptibility 
tensor in the laboratory frame Xg'^' ■ Then x(r;w) is the complex length of vector ^(r,a;) and 
'rt{uj) is its "director" , which are given by 

X^(r, w) = xl{r, uj) + x^(r, to) + xI{t^, ^) 

it = lt{uj) = ^(r, aj)/x(r, uj) (A9) 

We additionally assume that the atomic polarization is homogeneous along the atomic sample 
such that the "director" 'rt{u) is constant along the path in the phase integrals representation 
of the Green's function (lAGh . 

The local reference frame is linked with the laboratory frame via a rotational transformation 
characterized by Euler angles a, /3, 7. Then the susceptibility tensor components projected onto 
the plane orthogonal to the light ray direction and defined in the basis of circular polarizations. 
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are given by 



x+i+H-) 



(l + cos/3)2 1 (l-cos/3)^ _i . sm^P . 

-X+iv-) -I X-i\-) -I ;^ — Xov---J) 



X-i 



-1/ 



4 ' 4 



(l-cos/3)2^ +1, , , (l + cos/3) ^ 



2 

sin^ (3 



X-i 



X+i 



-1/ 



-e2'^sm2/3[x+i(...) + Xli(...)-2xS(...)] 
^e-2^^sin2/3[4}(...)+x:l(...)-2xS(...)] 



(AlO) 



where in the right side the tensor components Xq (•••) ^-re associated with the laboratory frame. 
In the basis of complex polarizations eo, e±i this tensor has a symmetric diagonal structure and 
its components are given by 



X^'(r,A) 



, no{T)_ 
" 2F + 1 



n(in),m' {m),m 



1 |(de*)„^p 
^ A + ir/2 



I ''ram' I 



1 



A + ir/2 Ac - A + S„„-(A) 



(All) 



where no(r) is the local density of atoms and the states |n) and \m') in the sum are unique for 
each |m) such that n = n{m) and m' = m'{m). 

The "crossed" term contribution in the differential cross section ()A2p is expressed similarly to 
the "ladder" term 



da 



12 



exp 



12 



i(k + k')(ri-r2) 



.{UJ12 - W2l) 

« ri2 



xXj,j(oo,r2,a;')aIfc' "^^ (a;i2-ki2V2)Xfc/(r2,ri,w^Ja;"''^"''' (a;-kvi)Xj^(ri, -oo,w) 



"^1' '^O"^' (^^21 -k2i vi)X^j(ri, r2, tJ2i)a^' (^-kv2)Xj^(r2, -oo,w) (A12) 
where for the reciprocal path the frequency 0^21 given by 

W2I = W - UJmqm2 + (k2i - k)v2 (Al3) 

is different from UJ12 but uj' for the scattering in the backward direction is the same. 

In the macroscopic and disordered atomic system the " crossed" term (1A12P contributes mainly 
in the narrow solid angle near backward direction (k' ~ — k). For any other directions this 
term makes negligible contribution. In higher orders of multiple scattering the "ladder" and 
interference "crossed" parts have the similar structure. All these term can be subsequently 
(moving from lower to higher orders) involved into the Monte Carlo simulation scheme. The 
cross section should be averaged with the atomic density matrix over the initial states mi, 
7712,. ■■ and over the atomic velocity distribution for vi, V2 ... This averaging is equivalent to the 
smoothed sum in the original expansion of the cross section given by Eq. (lAlh . 
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